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The energy loss of hard jets produced in the Deep-Inelastic scattering (DIS) off a large nucleus 
is considered in the collinear limit. In particular, the single gluon emission cross section due to 
multiple scattering in the medium is calculated. Calculations are carried out in the higher-twist 
scheme, which is extended to include contributions from multiple transverse scatterings on both 
the produced quark and the radiated gluon. The leading length enhanced parts of these power 
suppressed contributions are resummed. Various interferences between such diagrams lead to the 
Landau-Pomeranchuk-Migdal (LPM) effect. We resum the corrections from an arbitrary number 
of scatterings and isolate the leading contributions which are suppressed by one extra power of the 
hard scale Q 2 . All powers of the emitted gluon forward momentum fraction y are retained. We 
compare our results with the previous calculation of single scattering per emission in the higher- 
twist scheme as well as with multiple scattering resummations in other schemes. It is found that 
the leading (l/Q 2 ) contribution to the double differential gluon production cross section, in this 
approach, is equivalent to that obtained from the single scattering calculation once the transverse 
momentum of the final quark is integrated out. We comment on the generalization of this formalism 
to Monte-Carlo routines. 



PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 



I. INTRODUCTION 

The modification of hard jets in dense media [H-dj], 
and the use of this modification to measure the partonic 
structure of both confined or deconfined matter is now 
a considerably mature science Q. While the modifica- 
tion of full jet structure is only now being measured, 
the modification of the yield of leading hadrons, (a re- 
duction or suppression of the yield) has been measured 
by the HERMES experiment at DESY for Deep- Inelastic 
Scattering (DIS) on a large nucleus [f| and by both 
the STAR and PHENIX detectors [g| at the Relativis- 
ts Heavy Ion Collider (RHIC) for jets in heavy-ion col- 
lisions. On the theory side, there now exist at least four 
sophisticated jet modification schemes which describe the 
single particle yields within a pQCD based phenomenol- 
ogy. It has now become customary to refer to these as 
the Gyulassy-Levai-Vitev (GLV) scheme @, the Arnold- 
Moore- Yaffe scheme (AMY) Q, the Armesto-Salgado- 
Wiedemann scheme @ and the higher-twist (HT) ap- 
proach [IM1- Of these, the GLV, ASW and HT ap- 
proaches have also been applied to the case of DIS on 
a large nucleus using almost an identical formal- 

ism to that used in the case of hot matter. All these 
calculations obtained a much smaller value of the lead- 
ing transport coefficient in cold nuclear matter than for 
heavy-ion collisions. This represents one of the leading 
pieces of evidence for the formation of a new type of mat- 
ter in high energy heavy-ion collisions. 

The transport coefficient in question, q, defined as the 
transverse momentum squared gained by a hard parton 
per unit length, has been identified as the leading ef- 
fect of the medium on a hard parton Q, [l2| for the 
case of light flavor suppression. In the case of DIS on 
a nucleus, the q, assumed to scale with the nucleon den- 



sity, is a constant, as a function of time, over most of 
the space through which the parton traverses (ignoring 
edge effects in the Woods-Saxon distribution). The case 
of heavy-ion collisions requires more sophistication: q is 
usually assumed to scale with an intrinsic quantity such 
as the entropy density s oc T 3 (T is the local tempera- 
ture). In a dynamically expanding medium, T depends 
on both the location in space and time. Hence a com- 
plete calculation requires an underlying realistic model 
of the medium. Three of the formalisms (ASW, AMY 
and HT) have been extended to calculations where the 
evolution of the medium is treated with a 3-D relativistic 
fluid dynamical simulation [T7j . 

Current experiments have, however, progressed far be- 
yond a measurement of light flavor single particle inclu- 
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as well as triggered distributions of associated parti- 
cles. Theory calculations have also progressed, espe- 
cially in the case of two particle correlations. While the 
GLV " 



sive suppression and now include measurements of heav; 
flavor suppression and azimuthal anisotropy 
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and AMY [22j approaches have 
attempted to describe the away side associated yield us- 
ing a leading order partonic cross section, the HT has 
incorporated a Next-to-Leading Order (NLO) hard cross 
section (23|. Along with this, the HT approach has also 
remained the sole formalism to incorporate multi-hadron 
fragmentation functions in an effort to explain the near 
side associated yield [H| ■ 

In spite of these developments, the HT approach 
still suffers from the disadvantage that, unlike the 
ASW and the AMY schemes, it is a single scatter- 
ing scheme. In its original form [l0| . the HT ap- 
proach only contained one scattering and one emission 
in the medium. In subsequent efforts [2a |. the HT 
approach was extended to include multiple emissions 
by iterating the singlc-scattering-single-emission kernel 
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by means of a Dokshitzer-Gribov-Lipatov-Altarclli-Parisi 
(DGLAP) (2(| evolution equation. While this does in- 
clude the case where a hard parton or its radiated gluons 
scatter multiple times, these scatterings are included in- 
coherently, i.e., there is at most one scattering within the 
formation time of a single produced gluon. The applica- 
bility of this assumption can be quickly estimated: imag- 
ine a typical hig h energy DIS event (say at the HERMES 



experiment |27|) with a Q 2 ~ lOGcV 2 and a forward 
energy v ~ 30GeV. The produced quark is off-shell by 
at most Q 2 and this depreciates to a non-perturbatively 
small value in a series of iterative steps. Consider a step 
where the drop is AQ ~ 2GcV with a momentum frac- 
tion of 1/2, the formation length of this radiation is, 



2(AQY 



lfm. 



(1) 



Assuming that there is at most one scattering per nu- 
cleon, given that double scattering in a nucleon involves 
twist-4 matrix elements, and there is approximately one 
nucleon per fm of length, this particular emission will 
encounter one scattering during the formation of the ra- 
diated gluon. However, for smaller drops in Q 2 , smaller 
initial values of Q 2 and higher energies v, this simple 
approximation breaks down. The limit of at most one 
scattering per radiation cannot hold in the case of jets at 
RHIC either, as the density of scatterers is much larger 
than in cold nuclear matter. Such considerations require 
the extension of the HT scheme to include multiple co- 
herent scatterings per emitted gluon. It is the object of 
this paper to carry out this extension. 

This paper completes the efforts begun in Refs. fl2llT3l| 
by computing the single gluon emission cross section from 
multiple scattering of a hard quark in a dense medium. 
The prior efforts included the simpler problems of eval- 
uating the all twist contribution to the propagation of a 
hard quark produced in DIS on a large nucleus without 
gluon radiation and also the emission of electromagnetic 
radiation from this hard quark. Similar to the case of 
photon radiation at all twist in Ref. [HI], we will eval- 
uate the triple differential distribution to split the hard 
quark into a quark and a gluon separated by a trans- 
verse momentum l q j_ with the gluon carrying a large 
transverse momentum l± and a fraction y of the origi- 
nal quark light-cone momentum. In a real experiment, 
the size of the electromagnetic fine structure constant 
oiem ensures that multiple hard photon bremsstrahlung 
is suppressed, however given the size of a s at the typical 
scales involved, the produced quark will radiate gluons 
multiple times. In this effort, only the differential single 
gluon emission spectrum will be calculated. The inclu- 
sion of multiple emissions which involves an iteration of 
the single gluon emission kernel, followed by the eventual 
numerical calculation and comparison with experimental 
data will be presented in a separate effort. 

The paper is organized as follows: Sec. II will intro- 
duce the basic notation by evaluating the leading twist 
contribution to the differential collincar gluon radiation 



spectrum from a hard quark produced in DIS on a large 
nucleus. Sec. Ill will identify an inclusive class of di- 
agrams that yield leading contributions to the all-twist 
gluon emission spectrum. The collinear approximation 
will be invoked, integrations over the irrelevant momenta 
will be carried out followed by simplifications in the nu- 
merator structure. Sec. IV will carry out the various 
pole integrations and identify the source of the Landau- 
Pomcranchuk-Migdal effect. Sec. V will include discus- 
sions of the color factors and length enhancements of dif- 
ferent classes of all- twist diagrams. Sec. VI will evalu- 
ate the all-twist contributions in two different limits and 
compare with the calculations in the ASW and AMY 
formalisms. In Sec. VII we present our concluding dis- 



cussions. 



II. LEADING TWIST AND COLLINEAR 
GLUON RADIATION 

In this section we calculate the next-to-leading order 
correction to semi-inclusive DIS on a large nucleus with 
a quark in the final state. By next-to-leading order, we 
simply mean including one interaction term in the am- 
plitude and complex conjugate, which converts a single 
quark in to a quark and a gluon. As in the previous 
efforts [l2l [iH ] we will assume and refrain from an at- 
tempt to discuss factorization [28|, [2|| . We re-derive this 
straightforward result to familiarize the reader with our 
conventions and the choice of gauge. 

Consider the semi-inclusive process of DIS off a nucleus 
in the Breit frame where one quark with a transverse mo- 
mentum l q± and a bremsstrahlung gluon with transverse 
momentum lj_ are produced (see Fig.[T|), 

C(Li) + A(p) — ► C(L 2 ) + q(l q± ) + G(i_i_) + X. (2) 

In the above equation, L\ and L2 represent the momen- 
tum of the incoming and outgoing leptons. The incom- 
ing nucleus of atomic mass A is endowed with a mo- 
mentum Ap. In the final state, all high momentum 
hadrons (h\, h,2, • ■•) with momenta p\ , pi , . . . are detected 
and their momenta summed to obtain the quark jet mo- 
mentum (l q ) and gluon jet momentum (/) and X denotes 
that the process is semi-inclusive. 

Our choice of light-cone component notation for four 
vectors {p = [p + ,p~ ,p±\) is somewhat different from the 
regular notation, i.e., 



p + p 3 

p = — o — ; 



p 



p -p 



(3) 



In spite of the asymmetric definition, the formal struc- 
ture of all dot products will be identical to the standard 
notation. The kinematics is defined in the Breit frame as 
sketched in Fig. [T] In such a frame, the incoming virtual 
photon 7* and the nucleus have momentum four vectors 
q, P A given as, 



3 



q = L 2 -L 1 



2q- 



,0,0 



P A =A\p+, 0,0,0]. 



In this frame, the Bjorken momentum fraction of the 
incoming quark with respect to its parent nucleon which 
has the average momentum [p + , 0,0,0] is given as = 
Q 2 /2p + q~ . After the scattering with the virtual photon, 
the final state quark has a four momentum ~ [0, q~ , 0, 0] 
i.e., a large momentum in the negative z or negative light 
cone direction. This outgoing quark will radiate a gluon. 
The radiated gluon has a transverse momentum of l±_ and 
carries a fraction y of the forward momentum q~ of the 
quark originating in the hard scattering, i.e., 



(4) 
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FIG. 1: The Lorentz frame chosen for the process where 
a nucleon in a large nucleus is struck by a hard space-like 
photon leading to the production of an outgoing parton and 
a radiated gluon. 

The double differential cross section of the semi inclu- 
sive process with a final state quark (which will eventu- 
ally become a jet) with transverse momentum l q± and a 
final state gluon with transverse momentum l± 3> Aqc d 
may be expressed as 



El 2 cI(t 



dW^ v 



d 3 L 2 d 2 l q± d 2 l±dy 2ttsQ 4 a d 2 l q± d 2 l L dy' 



(5) 



where s = (p + L\) 2 is the total invariant mass of the 
lepton nucleon system. In the single photon exchange 
approximation, the leptonic part of the cross section is 
described by the particular form of the leptonic tensor 
denoted as given as, 



1 



= -Tr[£ l7M £ 2lv \. 



(6) 



In the notation used in this paper, \A;p) represents the 
spin averaged initial state of an incoming nucleus with A 
nucleons with a momentum p per nucleon. The general 
final hadronic or partonic state is defined as \X). As a 



result, the semi-inclusive hadronic tensor in the nuclear 
state | A) may be defined as 



W 



=Y,(^f5 i {q + P A -px) 



X 



x (A;p\J»(0)\X)(X\r(0)\A;p) 

I d 4 yo e^y°(A- lP \J^y )r(0)\A;p) 



(7) 



In the equation above, the sum (Ylx) runs over all possi- 
ble hadronic states and = Q q 4> q j >J "if' q is the hadronic 
electromagnetic current, where, Q q is the charge of a 
quark of flavor q in units of the positron charge e. It 
is understood that the quark operators are written in 
the interaction picture and two factors of the electromag- 
netic coupling constant a em have already been extracted 
and included in Eq. (|5|). The leptonic tensor will not be 
discussed further. The focus in the remaining shall lie 
exclusively on the hadronic tensor. This tensor will be 
expanded order by order in a partonic basis with one on- 
shell gluon in the final state and leading twist and maxi- 
mally length enhanced higher twist contributions will be 
isolated. 

The leading twist contribution is obtained by expand- 
ing the products of currents to next-to-leading order in a s 
to account for the radiated gluon. This contribution may 
be expressed diagrammatically as in Fig. Q] This repre- 
sents the process where a hard quark, produced from one 
nucleon in a deep-inelastic scattering on a nucleus, pro- 
ceeds to radiate a hard gluon and then exits the nucleus 
without further interaction. Diagrams where the hard 
gluon is emitted from the incoming quarks are suppressed 
in our choice of light cone gauge. In the following, we an- 
alyze this contribution in some detail. Indeed, there is 
no new information presented in the current section and 
the discussion of such contributions is now well estab- 
lished [26} . The approximations made in this section will 
form the basis for the analysis of single inclusive gluon 
radiation at all twist. The semi-inclusive hadronic tensor 
may be expressed as, 



TF, 



Ai iv 



ACfWg v 

I d 4 yo(p\^(y )r0 00 l^(0)\p) 



(8) 



In the above equation, Cp expresses the probability to 
find a nucleon state with momentum p inside a nucleus 
with A nucleons. In the collinear limit, the incoming 
parton is assumed to be endowed with very high forward 
momentum (po + = XqP + ,p$ — > 0) with negligible trans- 
verse momentum po ± <C Pq . Within the kinematics cho- 
sen, the incoming virtual photon also has no transverse 
momentum. As a result, the produced final state parton 
also has a vanishingly small transverse momentum (i.e., 
with a distribution S 2 (p±)). In this limit, the leading 
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spin projection of the pieces which represent the initial 
state and final state may be taken. The factors, 



7° + 7 3 



7 ~7 



(9) 



are used to obtain the spin projections along the leading 
momenta of the outgoing state and the incoming state. 
The coefficients of these projections are the two func- 
tions, 



F(y ) = A(p\ij(yo)^i>(0)\p) 



(10) 



dition, 



s[(po + q 
S[-Q 2 + 
1 s 

2p+q- 
5[x - x B 



-I) 2 } 

2 P +( q - -n-2 q n- 

zo(l - y) - x B (l - y) 
xl] 



2p+q-{l-y) 



2q-l+] 
I 2 

. J= 
2p+q-y 



(13) 



where, the collinear condition that po ± — > has been 
used to simplify the final equation. The new momentum 
fraction Xl has been introduced: 



and (in a notation where the superscript on the operator 
O 00 implies that the quark undergoes no scattering in the 
initial or final state) 



1 Q00 

2 

dH 



(11) 



x Tr 



(2tt) 4 

7" 



d 4 zd 4 



2 (po + q) 2 -ie 



, dH q d 4 p d 4 p' 
' (2tt) 4 (2tt) 4 (2tt) 4 

' z 7 Q y q 2nS(l 2 q ) 



x G a ,(027r5(/ 2 )(-^) 7 -^P 



(p'o + q) 2 + ie _ 



x giq-Va e -i{pa+q)-{ya-z) e ~U-{z-z') e -U q -(z-z') 
x e -iCPo+9)V fl 2 > 



Integrating over z and z' yields the two four-dimensional 
(5-functions: <5 4 (po + 5 — I — lq) and 5 4 (po — p' Q ). 

The reader will have noted that we have ignored var- 
ious projections such as those which arise from the (-L)- 
components of the 7 matrices, e.g., 



C^Tr 



7 



.7 



(12) 



This approximation may be justified in the high energy, 
collinear limit l\jy « Q 2 where such contributions are 
suppressed compared to those of Eq. ([5]). As pointed out, 
in this effort, we will often adjudicate the importance 
of different terms using scaling arguments inspired by 
SCET. The kinematic regime explored in this article is 
l± ~ XQ and both y, (1 — y) < 1. In this region, l\jy ~ 
X 2 Q 2 <?C Q 2 . Another interesting regime is the soft y 
limit where y ~ A; even in this region one may ignore the 
non-transverse projections of the virtual photon. In so 
doing, the focus of the remainder of this article has been 
limited to projections where the incoming virtual photon 
is transverse. 

The on-shcll 6- function over / is used to set / + = 
l\j2l~ . The other on-shcll ^-function, instills the con- 



/ 



XL 



2p+q y(l - y) 



P+Tf 



(14) 



where y has already been defined as the momentum frac- 
tion of the radiated photon (l~/q~) and t/ ~ 1/(X 2 Q) is 
the formation time of the radiated gluon. 

The factor G a p(l) in Eq. (fTTj) represents the radiated 
photon's spin sum. In this effort, the light cone gauge 
(A~ = 0) will be assumed, i.e., 
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n ■ I 



(15) 



where we have introduced the light cone vector n = 
[1,0,0,0] which yields I ■ n = l~ . Note that with this 
choice of gauge, the largest component of the vector po- 
tentials from the initial states may still be regarded as 
the (-f)-componcnts. 

Substituting the above simplifications in Eq. (jlll) . 
leads to the simplified form for the final state projection: 



QUU = 



x Tr 



d 4 po dl d 2 l± 
(2tt) 4 (2tt) 3 2Z- 

7" i + q 



2 2p+q~(x a - x B - x Do - ie) 

{7?7-(fa) - xb]p + - (-g± aP ) 
-p- 7-L ■ H ~ 7 7± ■ l —r=- 



+ 7 i + {q ~ l )i 



2l_ 



i + q 



x 2?r 



2p+q (x — xb - x Do + ie) 
S[x — xb - x L ] 



2p+q-(l-y) 



(16) 



The i5-function may be used to carry out the Pq = xop + 
integral; the absence of p$ and po ± from the integrands 
allows for these integrals to be carried out and con- 
strain the locations y^ ,yo x to the origin. Further sim- 
plifications may be carried out by noting that 7 ± anti- 
commutes with j±, while { 7 + ,7 - } = 21 (i.e., twice the 
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unit matrix in spinor space) and {7 ± ,7 ± } = 0. Replac- 
ing l~ = q~y, in Eq. (|16p , one obtains, 



III. GLUON RADIATION FROM MULTIPLE 
SCATTERING 







mi 



<%o + )<%oJ 



dyd 2 l±_ 



e ~i(x B +x L )p+y p+ 



(2p+q-)Hp+q-(l-y)xl 
2(x L p+ - 



y) 



2l 2 



a s f dydl\ 2 - 2y + y 



21+ 



2tt I 2 



!J 



(17) 



Rcintroduction of the final state projection 00 m 
Eq. ([5]), produces the well known and physically clear 
formula for the differential semi-inclusive hadronic ten- 
sor with single photon emission in the final state, 



dydl 2 ± = Cp^^QJq (x B +x L ){-g^) 



(18) 



In this section, we start by taking the part factorized 
from the hard cross section denoted as and expand 
it as a power series in the hard scale Q 2 . The leading 
term was denoted as O 00 in the preceding section and 
only included the process of vacuum radiation without 
any scattering. Here we consider the effect of multiple 
scattering on both the quark and the radiated gluon. As 
in prior attempts in this direction, we will not discuss 
the issue of factorization }28l - l3lj ]. The focus will be on 
deriving the physical leading effect of multiple scattering 
on the process of single gluon radiation. 

We start with the part denoted as G q '^ nm i.e., the 
term with N scatterings on either side and with the gluon 
radiated from location q,p on the quark amplitude and 
complex conjugate and scattering a total of n < N — q 
and m < N — p times. This is represented in Fig. [2] 
The Feynman integral for this diagram can be written 
down in position space, with the quark propagator in 
the amplitude (left hand side of the cut) from y • to y' i+1 
written as, 



In the equation above, f q {xs + represent the parton 
distribution function of a quark with flavor q and electric 
charge Q q in units of the electron charge e, in a nucleus 
with momentum Ap where the parton carries a fraction 
(xb + xl) of the momentum p,i.e., 



! q { X B + xl) 



A 



d y 

2ir 



-i(x B +x L )p + y 



x ^(p\i>(yoh + ^)\p)- 



(19) 



In Eq. (JTSJ, the factor P q ^ ql (y) = (2 - 2y + y 2 )/y is 
the quark-to-photon splitting function; it represents the 
probability that a quark will radiate a photon which will 
carry away a fraction y of its forward momentum. The 



g" 



projection 

In DIS experiments, one has an experimental handle on 
the Bjorken variable xb- In this article we focus on the 
region where xb ^1- As a result, given that l± ~ AQ, 
we may approximate /(is + ii) - /(ib). As the par- 
ton produced immediately after the hard scattering has a 
vanishingly small transverse momentum, the transverse 
momentum of the final produced quark is simply the neg- 
ative of the radiated photon's transverse momentum, i.e., 



-Z 



q± — !i. As a result, the differential hadronic tensor 
for the transverse momentum distribution of the final 



quark is given as 



~ cpnY,Q 2 q f?(x B )(-g?) 



dydl\d 2 l q± 



: ?±-Lp q ^(y)S 2 (l ± +l q± ). (20) 



2ttZ 



S(yi+i -Vi) = 



(2tt)4 



q'i+i 



(21) 



The quark propagators in the complex conjugate (right 
hand side of the cut) are identical with the signs of the 
ie and the i in the numerator reversed. We then shift 
the momenta as q' i+1 = q + J2]=oPi wnere p'o is the mo- 
mentum of the original incoming quark which is struck 
with the virtual photon and the other p'^s are the mo- 
menta of the gluons off which the final out going quark 
scatters. Thus we replace the variable being integrated 
as dq' i+1 —¥ dp'j. Similar replacements are made in the 
complex conjugate. In the the gluon propagators, the 
replacements are made starting from the cut end, e.g., 
the i th propagator momentum in the complex conjugate 
Qi = I — Y^JLi hi where I is the momentum of the gluon 
propagator that is cut and the fc,'s are the momenta be- 
ing brought in by the "soft" gluons off which the hard 
gluon scatters, this replaces dQi — > dki and similarly in 
the amplitude. 

Having made all these shifts in momentum we can 
write down the full expression for the Feynman integral 
corresponding to the diagram of Fig. [2j We quote the full 
expression first and then will describe the various terms 
within. The full expression for the final state outgoing 
quark and gluon which scatter a total of N times in the 
amplitude and complex conjugate, with the gluon being 
produced at q in the amplitude scattering n times and at 
p in the complex conjugate scattering m times, is given 
as, 
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N—m N—n 







NN 
q,p\n,m 



11 11 rf ^ ^4 j^m (o^ 1111^ c 



i=l fc=l 
p 



(2tt) 4 (2tt) 4 (27r)*(27r)*Alll 



x Tr 



i=0 



-*{(e;= ^)+ 4 

[(EUPi) + ?] 2 -« 



N—m- 

7«o n 

2— p 



J (2tt) 4 (2tt) 4 

[(EpoP 3 ) + 9-i+(ELy] 2 -« 



(22) 



x 27r5(/ 2 )27r5(/ 9 2 )G^ (m,n) /, 



n 

N — i 
-m— 

n 

i=0 



- V 7/3o 1 



#0+ ^ 



fc=JV— 
N-m-1 



[(Ef=ofO + ? - ' + (E?=i fc 0] 2 + ^ J fe = 9 1 [(Ef=oPi) + «] 2 + ie J + 9)2 + ie 

m N—n — 1 n 



i=0 3=1 



x e 



j=i 

JV— m— 1 



»(jv-m)-s E p* I -9+i- E fc i ) i y'(N- n) -\ l <i-[ E Pj -9+'- E fc j 



x ( A 



N—m 



-igf<A+ ( yi ) igt b *A+ (y' k ) J] A+ (0) J] fl/^ 6 ^ (CD 



k=N-n 



3=1 



i=l 



The gluon line G"°^° (this refers to combination of gluon propagators and scattering vertices extending from the 
amplitude side to the complex conjugate) is given as, 



-iG 



[*-E^*i 



-(-.9±)/3 iQl 2r 



(-.9±) Qm /3^ amb '" n 



2^ (-5i)/3'a' 



<LT a;b;_i 



[*-E?=^T 



.(23) 




AP 



y/ jl' 



y q+ i 



• p • 



r — 3f ifft 



y P+ i y D % y 



AP 



FIG. 2: DIS on a nucleus with a produced quark and an emitted gluon scattering N times. 

I 



The Feynman diagram corresponding to this equation is tegrals over coordinates and momenta. The coordinates 
Fig. [2j The first line of Eq. (|22|) contains a series of in- yi, where 1 < i < N — m, represent the locations where 
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the quark scatters in the amplitude; the coordinates y' k , 
where 1 < k < N — n, represent the scattering of the 
quark in the complex conjugate. The momenta Pi-\, are 
the incoming momenta into the quark line with po being 
the momentum of the incoming quark which is struck 
by the virtual photon. The remaining momenta from 
Pi to pN-m-i arc momenta of the gluons off which the 
quark scatters. The same is true for the p' k momenta in 
the complex conjugate. The coordinates Q (Cj), where 
l<i<m(l<j<n) represent the locations where the 
emitted gluon scatters in the amplitude (complex conju- 
gate). The momenta carried by each of these insertions is 
characterized by fcj (k'j) where 1 < i < m (1 < j < n) in 
the amplitude (complex conjugate). The two remaining 
momenta I and l q represent the final on-shell momenta 
of the exiting quark and gluon as in the previous section. 

The second line of Eq. (|22[) represents the quark prop- 
agators on the amplitude side. The third line contains 
the cut propagators and the full set of gluon propagators 

from «o to /3q in one expression G"°f° which is decom- 
posed in Eq. (|23|). The Greek superscripts represent the 
vector indices where the gluon is first emitted from the 
quark line and the Latin subscripts represent the color at 
this emission. The fourth line contains the quark prop- 
agators on the complex conjugate side of the cut. The 
fifth line contains the phase factors for all but the cut 
quark propagator. This last phase factor is separately 
expressed in the sixth line. The seventh line contains the 
expectation in the nucleus for the N + N gluon field op- 
erators in the amplitude and complex conjugate. Here 
\A) represents the nuclear state |32| . 

Eq. (f2U)) . contains all the gluon propagators for the 
emitted gluon which undergoes m scatterings in the 
amplitude and n scatterings in the complex conjugate, 



where, G„ 



represents the numerators of the gluon 



propagators after the emission both on the amplitude 
and complex conjugate side. In the rest of this section 
we will discuss various simplifications that arise due to 
our choice of kinematics for the radiated gluon: Zj_ ~ AQ 
and l~ = yq~ ~ Q where Q is the hard scale in the DIS 
process. 

Eqs. (|22I23|) arc not completely general and have a 
set of approximations already included. While calcu- 
lations are carried out in the negative light-cone gauge 
A~. = 0, there remain three other components of the 
gluon field in covariant products such as 7u i ^4-a < which 
arise with each quark-gluon vertex. As shown in a re- 
cent study [HI, as well as in Appendix A, in the case of 
the Brcit frame, the dominant component of the vector 
potential (in A~ = gauge), for gluons radiated from 
partons with collinear momenta in the (+)-direction, is 
the (+)-component which scales as A + ~ \ 2 Q compared 
to the (-L)-component which scales as A± ~ A 3 Q. 

The numerator of the terms in the second and 
fourth line of Eq. (|2"2"j) have the general form 
• ■ ■liii < li i lv i A'^ i t ai .... While this includes all four com- 
ponents of the quark momentum and the gauge field, 



the (— )-componcnt of the quark momenta is q i = q + 
J2jPj — 1 ^ Q an d dominates over the (_L) and (+)- 
components which scale as XQ and \ 2 Q respectively. 
This is also true of the gauge field as mentioned above. 
As a result, the leading term in each product of quark 
propagator and gluon insertion is . . . 7 + <7~7~ A + .t ai ... if 
the quark propagator is situated before the radiation of 
the outgoing gluon and is . . . 7 + (q~ — A + .t ai ... for 

quark propagators which appear after the radiated gluon. 
This is identical to the case of photon radiation from a 
multiply scattering quark produced in DIS [l3j . 

The major difference in this article from Ref. [li|, be- 
sides the scaling relations used, is the scattering encoun- 
tered by the radiated gluon. This is somewhat non-trivial 
in A~ = gauge. We begin by elucidating the various 
propagator and three gluon vertex structures. We will 
only consider three gluon vertices as four gluon vertices 
indicate a scattering of two gluon fields at the same loca- 
tion. As pointed out in Ref. [L2\ , such contributions are 
not length enhanced and are thus suppressed. From the 
effective theory perspective, such contributions are not 
contained in the leading Lagrangian J33[. 

The gluon propagator in A" =0 gauge for a gluon 
with four-momentum I contains the constant unit vector 
n M = (1, 0, 0, 0) and is given as (suppressing color indices) 



G» v (l) 



iG^ v 
I 2 + it 



-9^ 



l-n 



l 2 + it 



(24) 



All components where either of the Lorentz indices is (— ), 
are vanishing, i.e., 

-1 + 1 



G~ + = G+- = 



r 2 



o. 



G~ 



G J 



0. 



(25) 



The leading components (where the numerator has the 
largest power of A) are those where both Lorentz indices 
are transverse, i.e., 



l— 7?. ■ = l- 



l 2 + it 



(26) 



The meaning of the notation Zj_ = is one where 

the index, which may be referred to as p or ± M can only 
assume values fj, =_L /J = 1,2, and (Z±) M = l ±fi = —l± = 
-(Z ± )„ and (ffxV = = 9^ U = {gxY v = 

— S J -' ±j . The _L_L components are followed by the (+ _L) 
and the (++) components, where 



G++ = 



G +±i 

"W+i^~ l (l-)[l 2 + itY 
G++ _ 21+ 

WTTt ~ l i-[i 2 + it\' 



(27) 



The three gluon vertex contracted with the gauge field 
may also be simplified as, 

(28) 



bac 



Tf^{l,k)AZ{k)=gf 
x [g Pa (k - 21) ■ A c + A ca (l - 2k) p + A c[i {l + k) a ] 

* gf bac gp a (-2r)At + A ca r 9 + + A cf3 r 9 +. 
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When the three gluon vertex is used in combination with 
gluon propagators on either side, the dominant term in 
terms of counting powers of A from the contraction of 
Lorentz indices is the first term on the third line of 
Eq. (EH), i.e., 



■ G° 



-lb 



g fa i - 1 b i+1 c iA + 2l - 



(29) 



The approximation of replacing the numerators of all 
the gluon propagators with the transverse projector g± ftv 
is consistent for all except the first and the last propaga- 
tors which connect with the quark line. Here the Lorentz 
index that contracts with the 7-matrix may assume both 
a + or a 1 component yielding the same power count- 
ing for the complete expression. As a result, in Eq. (|23|) . 
the numerators of the first and last gluon propagators 
may contain either a (G +± ) MI/ or a g^ v '. The numera- 
tors of all other gluon propagators contain solely the sum 
over transverse polarizations to obtain the leading power 
counting. 

As a result of the above simplifications, the fields 
in Eq. (f2"2"j) have all been replaced by A + . The contract- 
ing 7^ matrices are replaced by simply The mo- 
mentum dependent part of the three gluon vertices have 
been replaced by simply 2l~ as indicated by Eq. (|29|) . 
There remain the simplifications on the numerators of the 
quark propagator, which essentially amount to replacing 
all the numerators with j + q~ expect for the propagators 
that surround the gluon emission vertex. The contrac- 
tion with either a (G + ) ,w or a g^ leads to three terms 
from the numerator of the three propagators connecting 
to the gluon emission vertex, 



V" 



7 + ? _ 7,.— ^ y -l + (q--n (30) 



l- 



+ t! 7± M (-flT)7 + (ff--r) 



\i=0 j = l 

A similar structure arises on the complex conjugate side 
as well. Invoking all these simplifications leads to a 
considerable simplification of the numerator structure of 
Eq. (EH). 

Following the stated rules of power counting, the de- 
nominator of the various propagators in Eqs. (|22l23p may 
be simplified and expressed in terms of the usual momen- 
tum fractions (ToL IT3I ] . For the first set of quark denomi- 
nators, up to the gluon emission vertex, we introduce the 
fractions, 



2p+q~ 



(31) 



While written without vector attributes, it should be un- 
derstood that p l ± for all i are two dimensional vectors 
transverse to the direction of propagation of the jet. Us- 
ing this, we may write the i th denominator as, 



+ /T~ 



Di< p = 2p + q 



< D 



XB 



(32) 



where Xi = p~l /p + ~ A 2 and is thus included in the de- 
nominator. Ignored arc terms which depend on p~ which 
introduce X 2 Q corrections to q~ ~ Q. For the denomi- 
nators beyond the emitted photon, besides the fractions 
x l D which are defined similar to Eq. (3TJ we introduce 
the momentum fractions which connect the final outgo- 
ing transverse momentum of the radiated gluon with the 
incoming transverse momentum which strikes the hard 
quark, 



l^l 2 -2y/±-SLoPj 
2p+q-y{\ - y) 



(33) 



Another set of transverse momentum fractions arise from 
the combination of the incoming transverse momenta 
striking the gluon line and those striking the quark line, 



2p+q- 



(34) 



Finally there are the momentum fractions that connect 
the transverse momentum of the final radiated gluon with 
the soft transverse momenta striking it, 



2l± ■ k\ |fcj_| 2 ■ J2j=i+i k± 



(35) 



The forward momentum fractions of these incoming glu- 
ons denoted as, 



pi 



(36) 



With these fractions, the denominator of the i th quark 
propagator after the gluon emission can be expressed as, 



A> P = 2pV(i-y) 



XB 




(37) 



The denominator of the i gluon propagator has a much 
simpler expression in terms of these momentum fractions, 



D 



9 _ 



-2p + q y 



E- 



u 



(38) 



To Substitute the above simplifications, requires that 
we separate the trace over spinor indices contracted with 
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the vector indices in the numerators of the gluon prop- 
agators in a numerator factor INT, and the remaining de- 
nominators of all propagators in a denominator factor 
D, which are coupled by the set of integrals over position 
and momentum, i.e., 

o™,n, m = J®(yMp)x-v, (39) 



where, D(y) represents all position integrals on the am- 
plitude and complex conjugate side of the matrix element 
and £>(p) represents all the integrals over momenta. 

The numerator factor, containing all the 7 matrices 
from the quark line and all metric tensors from the gluon 
line, keeping only the leading terms in A, simplifies as, 



Tr 



-(7 7+g ) p 7 { 7+q 7m - 



h E7=i *i 



-7+(T-n + 7 2(i>i 7x M (- ff r)7 + (g--r) 



-r^i-gThl (5>Y + E fc i-^| } {T7 + 0r i-)} N ~ m - '-Vt+Ot r) 
=0 3=1 



{l-l + {q- I- )} N - n - 9 " 1 7" ^ 7 + (<T - r) 7p 



(40) 



+ 7 + (?--I-M(-ffiU Ei/i 



1 

■J* 



Ki=0 



£yi +E fc 'i-^ ) l P A-9t) P vl + 1-\ (7-7V) 9 

:0 j = l 



(2g-)"7 _ < 7 + T7„- 



- 7 + ( 9 --r) +71^1 7i M (-9r)7 + (?--n 



p 



u-<?r hi ( E^ + E fc i - 1 -l ) ) n} N - m - p - 1 2( q - n 
=0 j=i 



x {2( ? - - r)}^-"-'-^- |7 + (<r - r) 7p - 



-7 9 



(41) 



E/l +E fc 'i - ^ ) 7l(-3±)p,7V J> (2<T) 5 

:0 J=l 



The entire set of transverse components of the metric tensor {—g± tiv ) that arise from the propagators and vertices of 
the radiated gluon line lead to the contraction between the v component of the gluon emission vertex in the amplitude 
with the same component in the complex conjugate. The overall minus sign is due to the fact that there is one more 
propagator in the gluon line than vertex i.e., the cut propagator. 

The equation above for Ji does not contain any (+) or (— ) components of the incoming momenta: Pi, p[, k%, k[. 
The only non-transverse components are the large (— ) components of the outgoing quark (q~) and radiated gluon 
This allows for the integrations over all (+) and (— ) components of the incoming momenta to be carried out 
solely by considering the denominators of the various propagators. These set of integrations are carried out in the 
next section. 

The remaining expression for the numerator of 0^^ n m can be further simplified to obtain a more suggestive form. 
Isolating the coefficient Cl of (l± — E"=i ' C-L — Ej=i yields four terms, which may be summed to obtain, 



Cr 



(2<r) 



2N-m-n+l 



(2?-) 



y 

2N -rn-n+l 



■0--v) 



2N —m—n—p—q— 1 



"■{l-yf + 2(1 -y)+2(l-y) + 2y 



■(i-y) 



2N — m—n—p — q— 1 



y 

-(2-2y + y 2 ) 

y 



(2q 



-\2N-m-n+l 



"(1-2/) 



2N —m — n — p—q— 1 



(42) 



my)- 
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The coefficients of y J% =1 Px ■ (l± - E"=i ^1) . U± - E™ 1 - V E*=i and y 2 ELi Pi 1 E*=i P'i arc similar 
and may be included with the contributions containing Cl to obtain the general numerator as 

(n„-\2N-m-n+l ( m P \ ( n 1 

N = [zen (1 _ y) 2N- m - n - P - q -i 2P(y) i r± _ j^ki - yj^pi )■[ i± £ fc'i - »E^x I ■ ( 43 ) 

^ V i=l i=l / y j=l 3=1 

Ignoring the factors of y, (1 — y) and g~ on the right hand side, the remaining expression for the numerator 
has a very simple interpretation. While P{y) is simply the regular vacuum splitting function, the two terms, 

(j± — EI=i ^±~U Ei=i Pl) an( i ('-L — E™=i — 2/ Ef=i arc the transverse momentum of the radiated gluon 
immediately after the radiative vertex in the amplitude and complex conjugate respectively. In the next section we 
will simplify the denominator by carrying out a series of integrations over the light-cone components of the incoming 
gluon momenta. The numerator factor in Eq. (|43[) does not contain these components and thus plays no role in these 
integrations. 



IV. THE POLES STRUCTURE FOR SINGLE GLUON EMISSION FROM N SCATTERING 

In this section, the integrals over the light cone components of the incoming gluon momenta will be carried out. As 
mentioned in the preceding section, for this part of the integration, we only need to consider the denominators of all 
the propagators. The phase factors for the last scatterings on the quark at locations y/N—m) may be simplified, by 
the introduction of an (N — m) th momentum, using the four 5-function as 



// /N—m—l \ m 

d i p(N-m)S i I lq - I £ Pi) - q + l- ^k 



P(N-m) ■ (44) 



Using this the ± components of this four dimensional (^-function and the various momentum fractions, the argument 
of the 5(lq) function may be simplified as, 

N-n ro \ . § (l+ — 1 -£l\ / N — n 



r / N-n m \ . § [1+ — 2i / N-n ro 

J \ i=o »=i / J *h \ l=0 i= i 

* ' ,6 ( N f Xi - ^ +( t } ' - *b - A N ~ m) + 5> + ^) ■ (4 

2 P +q-(l-y) (1-y) J (l-y)J [ 

We also perform the integration over l + using the on-shell delta function, 

dH n _, l2 , /• dl-d 2 l ± 1 



We make the substitution, Z~ = <7~y and express all further expressions in terms of the radiated momentum fraction 
y. This simplifies the last two phase factors in the expression for 0^^ n m at locations y(N-m) and V(n— «)• Including 
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only the denominators, (^-functions (associated factors of 2n and i) and phase factors from Eq. [22j we obtain, 



N—m N-n 



NN 
q : p\n : m 



n-pr 4 4 , dx l -id 3 p i -i dx'^d^p'^ dyd 2 l ± dx {N _ m) (fp (N _ m) " " dzjd^kj dz^Pfy 



i=ij=i 



(JV-m) m 

i=o j=i / i=o I 2g 



AT-m-l 



n 



J I 2,- (1 - ,i [(Ej =0 a, - ^f^) .r„ 4 + (V»' ; . 



' N-r 



(p+) 2 2 ? -(l-y) 

n 



£) T v-^T J , (N—m) 

x j — ^ xj% Xj^ 



i=0 



(1-!/) 



j'=i 



_ n _! I 2g-(l - y) ( Eto I' - "' D '" (X ' ?)i 



^O^Z 1-y 



xb -x'i + [ Er=i 



n 



fc= g I 



2a- 



n 



n 



fc]M-23-y[Er=i^-^- 



=1 -2g y 

N—m m N — n—1 n 

J~J e -*P<-(lK-tf(jV-n)) J^J e - ifc J-(0-aJjv-„)) TT^ e *Pi-(j/<-I/(N-„)) J^J g^-tCj-yJjv-™)^ 



i=0 



J = l 



i=0 



i=i 



(47) 



The extra factor of (l/2p + ) 2 is due to the fact that 
there are two extra propagators than integrals over (+)- 
components of momentum. The gluon propagators have 
an over all negative sign compared to the quark propa- 
gators. This is due to convention of calculating the mo- 
mentum in the i th gluon propagator in the amplitude as, 



-2p+q y 



2^j=i Z 3 y U 



(48) 



The reader will note that besides the phase factors, 
the ( — )-componcnts of the incoming gluon momenta have 
been neglected from the integrand. As a result, we can 
integrate over the (— )-components of all the incoming 
gluon momenta isolating the entire process on the nega- 
tive light cone, e.g., 



dPi 
2ir 



-"i-Pi (Vi 



-v'tn- 



y {n- 



))• (49) 



y'tN 

to the y' N 



Like wise integrating over the set of momenta p[ , k, t and 
k[~ yields the sets of 5-functions [S(y'f— y'ov-m))> 

)j HC't ~ y'(N-m)) 1 tnat constrains every vertex 
surface. The <5-function that arises from 
the integration over p' , equates, y'%_ m and thus all 
other (+)-locations to the origin (j/'q = 0). This removes 
all negative light-cone components of momentum from 
the derivation. It should be pointed out that, in the 
higher twist formalism, these components are responsible 
for elastic energy loss [HI, HH . By integrating out these 
components wc specialize to the case of only radiative 
energy loss. In reality, as is obvious from the formulation 
of the Feynman integral of Eq. (|2"2"|) , elastic and radiative 
loss are coupled, with each contributing to the magnitude 
of the other. We leave the coupled calculation of elastic 
and radiative loss to a future effort. 



The integrations over the momentum fractions Xi, x^, 
Zi and z[ requires a certain amount of care as regards the 
order in which these have to be performed in order to 
have the simplest pole structure. We focus only on the 
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remaining light-cone components of the phase factor, The remaining delta function over the (-f)-components of 

the momentum may be used to set the last momentum 

N-m m <V„„+;„v, .(— ^+ 

r + = jj e -^p (»« -y (w _„)) jj , 

i=0 3=1 
N-n-l 



ix iP + (yr- v %_ n) ) n fi -<. JP +(C i --v'^_ B) ) fraction *(JV-m)(= Pfr_ m) /p + )- Similar to the case of 

photon production and multiple scattering in Refs. |13| 
and Qjl, we obtain, 



i=o j=l 

(50) 

I 



x (N-m) 



X B 



(N-m) 



E 

^=1 



N-m 

E 

i=0 



AT-m-1 

E ^ 

i=0 



(51) 



This converts the light-cone phase factor to, 



N — m m N — n — 1 n 

Jx^p+iy'- -y'~ N n) ) TT iz' j p + (C'j -J/'(^_„) ) 



n 



x exp 



e -«iP + (?/i -«(jv_ m )) JJ g-^jP+tCj -y (JV - m) ) J~J e «b + (?/'i -j/' (n _„)) TJ , 

3=1 



i=0 j=l i=0 3=1 

3/ r) ""("" 2/< 



N-m 



W [y (N - m) - V (N-n) )\XB + 



(N-m) 



1 - y f-' T 



fc=i 



(52) 



The remaining integrations over momentum fractions involves contour integrations around the various poles of the 
2N + 2 propagators. The simplest means to carry this out is to start by integrating the penultimate momentum 
fractions on either side of the cut quark line, X(^-m-i) an d a;'(jv-n-i) an d proceed towards the lower momentum 
fractions. In this way, the fraction being integrated over appears in at most one propagator. The integration over 
X(N-m-i) requires the contour to be closed at -Moo and, as a result, the positions V( N _ m \ arid y7 N _ m _ 1 - > need to be 
ordered with y7 N _ m \ > y~ N _ m _iy The isolated residue of this integration yields, 



/ 



dx 



(N-m-l) 



-lX(jv_ m _i) 



-1) »(IV-n); 



\Lj=0 X 3 - (1-y) 



ie [y(N-m) ~ y [N-m-l) ) eX P 




(N-m-l) 
XB - XI 



l P [y(N-m-l) V(N-n) 



N-m-l 



(N-m-l) 
Xb + x\ 




(53) 



Similarly, the integration over x^ N _ n _^ requires the contour to be closed at — ioo and thus orders y'^jy-n) -> v' (N-n—i) ■ 

Integrating over single poles such as the above may be continued until x p +± in the complex conjugate (c.c.) and 
x' q+1 in the amplitude. At this point, there are two double poles in the c.c. in x p and z\, and two double poles in 
the amplitude in x' q and z[ . Eq. 1471 after all the contour integrations up to x p and x' q have been performed may be 
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expressed as (with integral signs implied), 



33 



NN 



-- [ [ <ryi<Ty k 



i=i k=i 



(2tt) 2 (2tt) 2 (2tt) 3 j/ (2tt) 2 



i dz%d k dz jd k \ 
j ~{2nf (2tt) 3 



(iV-m) 



x d% ±< 5 2 h, ± + ix - E - E fc i M I 



X < 



3=1 



2tt 






2q- 


i J2 x i~ x b ) ~ x b - ie 






.v=° / 





2?"(l-y) 
l 



E - 

V 3=0 



_ ~-P 



E • 



i 



(p+) 2 V2g-(i-y) 



N—m—p— 1 



l 



2<r(l-i/) V2g-(l-y) 



1 



i-v 



N—n—q—l 



dxi 



2ir 


(jtx'i-x'n) 




2q- 


- x B + it 









2?-(i-y) 



9 
1=0 



X 



II — 

i=i -2g-?/ 

n 

E=l 
—7) 

n 



Zk - ^tT~ + le 



n — 

J_1 -2g-» 



1=3 



N—m—1 m TV— n— 1 n 

e vi-(»i-v'x v_n) ) TT eifci-CCi-s/'f-"') JJ ft -»p'x-(w'i.-»'x v ~ n) ) TT fi-^'i-CCWf"* ) 



i=0 
N—m 

x | e 

z=p+l 



3=1 



i=0 



n 



n 

P 9 
TJe^P+^y,— Vfi) JJ e iP + x'j{y'J-y'- +1 ) 



3=9+1 

m n . _f_ _ 

JJ e -ip + ^(CT-y- +l ) JJ e ip + ^(C'7-y'- +1 ) xe " lp «W) 



1=0 



3=0 

-I J^^Sr)] 



3=1 



x e 



(54) 



r 



The px , J/x aie two dimensional vectors with a Euclidean 
metric. The new momentum fractions, Sx z L and Sx' J L are 
defined as in Refs. 1131. 



2 2—1 

X L ~ 



2p+q y(l - j/) 



-2Zx • p\ 



8x\ 



2p+q-(l-y) 1-y 



(55) 



The momentum fraction on the amplitude side, Sx L , has 
an analogous definition and like 5x l L depends only on 



the transverse momenta of the gluons which strike the 
propagating quark: 



T I'J _ 13-1 
x L x L 



-2yl±-EL p' l ± + 2yl±-EUp' 



2p+q y(l - y) 



-2l ± ■ p 



t.i 
1 



Sx'l 



2p+q-(l-y) 1-y 



(56) 



The first set of non-trivial integrations are those over 
x p and x' q , which involve two poles. It should be pointed 
out that this is not a true double pole and are indeed 
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two separate single poles: one may simply take the sum over x p , we obtain, 
of the residues at the two poles. Isolating the integration 



J 



-ix p p+(y p -y (p+1) ) 



— I 



2?r 



E^i - x d —xb — it 

i=0 



P 

E X i 

i=0 



-^(x£ + {| o ^^}-{| /i+ W})(^-^ +1 )) 



(57) 



There is a completely analogous expression for the integral over x' q in the amplitude which involves the integral over 
two poles. The remaining integrations over the x$ for i < p and over x'^ with j < q are rather trivial and consist of 
simple contour integrations as performed in Refs. [HI, [35j]. When the set of propagators and associated phase factors 
in the equation above is combined with the propagators and phase factors in the gluon propagators there exist two 
poles for the momentum fraction z\ in the complex conjugate and z[ in the amplitude. Isolating the integrand of the 
z\ integration, we obtain, 



/ 



dz-\ 



i—1 



E. 

! = 1 



y> V x d + ( x t ) j 

2^ i-j, 

i=0 



-ip+*i(d -l/„ +I ) 



(*?)< 



U=l 

—if 



yx' D + (x!py 



i=0 



X < 



E *»+E 

t=2 i = l 
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*p + - E «i+E 



U=0 1 



+ %g+i - Ci >e 



(- E *-§ ^2i+<4+{ E ^?L_}) ( Cl -»- +1 ) ' 



E 



2 = 1 1.1=0 



E 



1-3/ 



* + e^-e ^ (cr-»p + i) 



i=2 i=l 



0(Ci - ) {* 



»p + < m (y P -»p+i) _ j e -*p < m (Ci -fp+i) 



(58) 



Where the aptly named momentum fraction x v ™ is expressed as 



yx), + (x$)i 



^2 y d 



E 



i± - E fc I - y E Pj 

i=l j=0 



2p+g j/(l - y) 



(59) 
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Similarly, the integral on the amplitude side is given as 

-( fc E<-^)(cW- +1 ) 



■ ^(CV - v7+i) - if (i/'^+x - c7)e ip+a: '* n «W g - +1 ) 

where the momentum fraction on the amplitude side is given as 



(60) 



ii-En-i/EfT 

i=l fe=0 



2p+g ?/(l - y) 



(61) 



The remaining integrations over the variables Z2 to z m , z' a to z' n , x p -i to xq and to can now be performed 
to complete all contour integrations. The resulting expression for T)^ n m is now completely ordered in the negative 
light-cone direction, i.e., all the y~,y'~ , C~,C~ locations are strongly ordered: 



N-raN-n rf2 i-l J2 l k-1 , ,2/ v JN-rn) m n 



0^,*= I I <*W 



Vk ' 



i=l fe=l 



(2tt) 2 (2tt) 2 (27r) 3 y (2tt) 2 



j2p ^2jL/fe 



j (2tt) 3 (2tt) 3 9± 



i=i j=i 



/ (N-rn) m \ , 



p+i 



l l 



2«-(l-i/)xr2p+- \2q-(l-y)J 2q-(l-y) 



N—m—p—l 



N—n—q—1 



1 1 



2q-(l-y)J 2p+ x'f2q-{l-y)\2q- J \2q~ J \2q 



9+1 



x - i/p ) {e-* + *r»,- _ e -<p + xrcr } _ 0(Ci - _ y - +i)e -ipM m % +I _ _ Cr)e -^ + *rcr 

x >(CT " 2/7) {e* +<e '^*'« - e^TCT } _ gg- _ ^Je^'*"^ - - C7K P+ *' rC ' 

p N — m 



fa l+ipi-yi 



4=1 



i=p+l 



x sec > c-_x > > on e - ip+cr ^ +ifei< ^(c; > Cn-i > > co n e^v^-^i^ 

i=i j=i 

9 JV— n + , 



-zp ^-J, • 



(62) 



In the equation above, the exponentials containing the dot products of transverse positions and momenta have 
changed. We have defined a new transverse momentum, 



-N-r 
P'± 



N- 



m m 



N-n—l n 
i=0 j=l i=0 j = l 



(63) 



Since there is no integration over y'\, we may shift the p'\ integration as d 2 p'°j_ — > d 2 p'^_ ". This allows the nesting 
of the various interactions, where transverse momenta from the amplitude and complex conjugate may be combined 
in pairs. This has to be done in a way so as to obtain the maximal enhancement in length and color and will be 
carried out in the next section. ^ 



V. COLOR FACTORS AND LENGTH 
ENHANCEMENT OF MATRIX ELEMENTS 

With the completion of the contour integrals over both 
light cone components of the incoming momenta and the 



simplification of the numerator Dirac structure, there re- 
mains the sum over all interactions and the simplification 
of the color factors. Unlike the case of photon radiation, 
the color factor in this case is very diagram dependent 
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and depends strongly on where the radiated gluons at- 
tach to the quark line as well as the number of cross 
connections between the scattering of the gluon and that 
of the final quark. The determination of the color fac- 
tor has to be carried out in tandem with the estimation 
of the length enhancement achieved by each diagram. 
As a result, in contrast to the case of photon radiation 
the first step is to decompose the expectation of the 2N 
gluon operators in the nucleon state. This is followed by 
an evaluation of its color factor and length enhancement. 
Only after this can the sum over emission points of the 
radiated gluon be carried out. 



The decomposition of the nuclear state as indicated 
by Fig. |31 may be carried out as in Refs. [H, 0: The 
nucleus is approximated as a weakly interacting homo- 
geneous gas of nucleons. Such an approximation is only 
sensible at very high energy, where, due to time dilation, 
the nucleons appear to travel in straight lines almost in- 
dependent of each other over the interval of the interac- 
tion of the hard probe. In a sense, all forms of correlators 
between nucleons (spin, momentum, etc.) are assumed 
to be rather suppressed. As a result, the expectation of 
the 27V gluon and 2 quark operators in the nuclear state 
may be decomposed as 



A. Color factor and length enhancement for the 
fully symmetric diagram 

Given the different color factors for different diagrams 
we need an organization scheme. We begin with the 
completely symmetric case where m = n < N and 
p = q < N . In this case we may nest the various in- 
teractions as indicated by the diagram in Fig. [3] Each 
blob represents one nucleon off which the quark or gluon 
may scatter. As in the case of transverse broadening [l2| 
and photon production [l3||, we will ignore the expecta- 
tion of more than two operators per nucleon state. We 
will not explicitly consider the expectation of two succes- 
sive gluon operators on the same side of the cut line in 
the same nucleon state. These do not induce any change 
in the transverse momentum distribution and constitute 
unitarity corrections. They are included in the final re- 
sult by a normalization of the resulting transverse mo- 
mentum distribution 11211. 




FIG. 3: DIS on a nucleus with a produced quark and an emit- 
ted gluon scattering TV times symmetrically in the amplitude 
and complex conjugate. 



2N 

(A-p^iy-^x^+^llA+iy^A-p} 

2N 

= ac* ^m y -, v ±)7 + m n a+ ( yi )\ P i) 

i=l 

+ C P l P2 (Pil^(2/~ 5 ^)7 + ^(0)bi) 

27V 

x ( P2 |n<(^)b2> + ..., (64) 

i=l 

where, the factor represents the probability to find 
a nucleon in the vicinity of the location y, which is a 
number of order unity (it is the probability that one of 
A nucleons distributed in a volume of size cA will be 
found within a nucleon size sphere centered at y). The 
remaining coefficients represent the weak position 

correlations between different nucleons. The overall fac- 
tor of A arises from the determination of the origin (the 
location in the equation above) in the nucleus, which 
may be situated on any of the A nucleons. Solely for the 
current discussion, we reintroduce the quark operators 
^{v~iV.l) an d V'(O) in the above equation. 

It is clear from the above decomposition that the 
largest contribution arises from the term where the ex- 
pectation of each partonic operator is evaluated in sepa- 
rate nucleon states as the j/i integrations may be carried 
out over the nuclear volume. As a nucleon is a color 
singlet, any combination of quark or gluon field strength 
insertions in a nucleon state must itself be restricted to 
a color singlet combination. As a result, the expectation 
of single partonic operators in nucleon states is vanish- 
ing. The first (and hence largest) non-zero contribution 
emanates from the terms where the quark operators in 
the singlet color combination are evaluated in a nucleon 
state and the 27V gluons are divided into pairs of singlet 
combinations, with each singlet pair evaluated in a sepa- 
rate nucleon state. This reduces the discussion to TV + 1 
nucleons, one in which the original quark is struck by 
the virtual photon and the remaining N off which the 
struck quark scatters on its way out of the nucleus. As 
has been demonstrated in Rcf. [36| . under the assump- 
tion that there are no correlations between the nucleons, 
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the combinatorial factor reduces to, 



C A = AC A ( - 

K - y p 1 ,...,p N+1 ^ w pi I 



N 



(65) 



where, Cp\ varies with the nuclear density profile p and is 
equal to unity within the nuclear radius for a hard sphere 
distribution. 

Further simplifications arise in the evaluation of gluon 
pairs in a singlet combination in the nucleon states by 
carrying out the y± integrations. The basic object under 
consideration is (ignoring the longitudinal positions and 
color indices on the vector potentials) 



dY±d 2 y' 3 ± (p\A + (y\)A + (y' J ± )\p) 



x e 

= (2tt) 2 (5 2 (p! -p'i) J d 2 y^e- ixip+(v t- y 'i ) 

x e^-yHp\A + (yj2)A + (-yj2)\p), (66) 

where, y± is the transverse gap between the two gluon 
insertions and p± = (jp\ + p J ± )/2. The fractions re, and 
x\ represent any of the possible types of momentum frac- 
tions xd,xt,xl and any combination thereof which ap- 
pear in the expression for 0p^ m m . The physics of the 
above equation is essentially the transverse translation 
symmetry of the two gluon correlator in a very large nu- 
cleus. 

One will note that the two dimensional delta function 
over the transverse momenta has removed an integration 
over the transverse area of the nucleus thus reducing the 
overall A enhancement that may be obtained. This is 
then used to equate the transverse momenta emanating 
from the two gluon insertions in the amplitude and com- 
plex conjugate amplitude. This also simplifies the lon- 
gitudinal phase factors which now depends solely on the 
longitudinal positions of the two gluon insertions. The 
assumption of short distance color correlation length (an 



assumption valid both in the deconfined plasma as well 
as in color confined cold nuclear matter) leads to the con- 
straint that the longitudinal positions of the gluon oper- 
ators from the amplitude and complex conjugate, whose 
expectation values are sought in the same nucleon state, 
cannot be more further apart than the color correlation 
length, i.e., 



dy-dy' (p\A+(y-,y x /2)A+(y' ,-y±/2)\p) 
dy-y- (p\A+ (y- , j?±/2)A+ (y~ , -y ± /2) \p) (67) 



In other words the two locations y~~ and y' have been 
restricted to be in close proximity of each other and each 
single nucleon expectation of two gluon operators may 
introduce at most one factor of length enhancement from 
the remaining y~ integration. 

The requirement that the two gluon operators be in the 
color singlet state, may be enforced by averaging over the 
colors of the two gluon operators, 



(p\AtA+\ P ) 



* IP) = jNf^l) m+aA+M - (68) 

A similar constraint on the quark antiquark cxpection, 
yields, 



(p\il>i'Y + il>j\p) = jj-{p\^i + Mv)- 



(69) 



The expectation of the quark operators as well as the 



associated phase factor exp 



-ix B p + y 







PI -VI 



and 



combinatorial factors AC^ may now be extracted and 
combined to form the nuclear structure function of the 
struck quark. Given these simplifications, the full expres- 
sion for the multiple scattering of the struck quark off N 
gluons in the symmetric pattern of Fig. [31 Qp pm,m> ma y 
be written as, 
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JV- 
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d 2 p\ dyd 2 l±_ 
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2e 
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%-i) e 



6 a * a * (p\A°* A°> \p) 



II {%r-yr-i)e 



i=p+l 
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AT- 



IV? -1 

xc>c-,>->g)n^ ^:f w 

=1 c 
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Jt a H b " ' ' t a ' t<t b '«Y[t 

i—1 i— p+1 i—N — m 



i—p 



i=l 



(70) 



where, (p|A 0i A 0i |p) is short hand for 

d5y i e- <J,+x o^. : " (71) 

x ( P \A°>(y- +6yr/2,y i L /2)A°'(yr - 5yr /2,-y\/2)\p). 

The overall color factor for the Feynman rule above may 
be evaluated by calculating the trace over the color matri- 
ces and dividing by the factors obtained from averaging 
over colors, which yields, 



C 



p,p;m,m 



(Cf) 



(JV-m+l) 



N C (N 2 - 1) N 



V ^ ) 



(JV-m+l) 



jV c (jV c 2 - l) w 

i / 1 \ iV-t 



jV. 



1 



2jV c 



1 

2iVc 



A, 



iV^-1 



.(72) 



Under the assumption that the parton propagates 
through a long extended medium, matrix elements are 
enhanced by the length traversed through the medium. 
Without any constraints except those of Eq. (|67|) which 
connect the scattering locations on the amplitude with 
those on the complex conjugate, the maximal length en- 
hancement factor is 




= L 



N 



(73) 



One can isolate the leading length enhanced contribu- 
tion for a given configuration by analyzing the various 
^-function constraints that limit the length integrations. 
For the diagram of Fig. [3] corresponding to Eq. ([70"|) , we 



note that all the scattering points on the quark and the 
gluon are ordered in the (— )-light-cone coordinate i.e., 
Cm > • ■ • > Cf' VN-m >•••>%• In addition the gluon 
scattering points for the first term in square brackets are 
forced to be ahead of y~ on the quark, i.e., £-f > y~ . 
For the second term in square brackets Cf > Vp+ii while 
for the third term y p+ i > Ci > Up ■ 

The largest length factor for each of the terms may be 
estimated by performing all integrations over the entire 
length with a unit integrand. For a single stream of con- 
strained integrals, this procedure yields the well known 
result, 



L t-L 

dy Q / dyi 

Jya 



dy. 



N 



VN+l 



N 



(74) 



For the case of two separate streams joined at a point 
y p as in the case the first term in the square brackets of 
Eq. (|70[) . this procedure yields the length enhancement 
factor, 



-AT, JV 

J p.p;m,m;l 



(N-p)\(j L dy) N 
(N — m — p)\(m)\N\ ' 



(75) 



To make it easier to analyze these factors we make the 
replacement, m q = N — m — p (m q is the number of 
times the quark scatters after radiating the gluon) , which 
yields, 



-JV.JV 



(ni + m q )\ [J dy 



N 



m q \m\N\ 



(76) 



Similarly, the length enhancement factor for the second 
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term in square brackets is given as, 



L 



N,N 

p,p;m,m;2 



(N-p-l)[[j L dy) 
(N — m — p — l)!(m)!AM 

{m + m q - 1)! (J dyj 
(m, - l)!(m)!iV! 



(77) 



Finally the length enhancement factor for the third factor 
with the square brackets of Eq. (|70|) may be obtained by 
noting that is locked between y~ and 2/^+1 > thus this 
is similar to the the first length enhancement factor but 
with the replacements m — > m — 1, p — > p + 1, i.e., 



L 



N,N 

p,p;m,m;3 



(N-p-l)\[f L dy) 
(N - m-p)\(m - l)\Nl 

(m + m q - 1)! (Jo dyj 
(m 9 )!(m - 1)!AT! 



(78) 



Inspection of these terms suggests that the largest con- 
tributions arise from the case where m + m q — > N, i.e., 
most of the scattering happens after the gluon has been 
radiated. We can now complete the analysis of this term 
by expanding the transverse dependence of the hard part. 
In this case, the hard part consists of everything except 
the matrix elements and the associated phase factors 
which contain the transverse momentum controlled by 
the matrix elements. The first step is to institute the 
condition that lj_ 3> kj_. This simplifies the momen- 
tum fractions which appear in the phase factors to 
Xl = (2p + q ~y(l — y))- One can also expand the sum 
of transverse momenta in the denominator as, 



term above and the transverse momentum dependent 5- 
function yields, 



n 



1 o 2 o 2 



±2 d 2 k\ d 2 v \ 



1 M£!±M£^ 

/2 + 74 



2l± • E fcj 
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(80) 
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2ml 1 • V; 



N-l 
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<-/± 



4m ( V 2 



JV-l 



<5 2 i 
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n i^i 2 i^ 



"II 2 - 



The factors of incoming transverse momentum, both 
on the gluon line (|&jj 2 ) as well as those on the quark 
line (|f^_| 2 ) may be combined with the position depen- 
dent matrix elements to yield the Lorentz force correla- 
tor obtained in Ref. [l2[ (suppressing color indices and 
longitudinal momentum and position arguments): 



1 



1 
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EI 



= 1 K _L 



j/E?=iPj 



(ELifcl + ^ELiPi 



/4 



(79) 



Since contributions are dominated by smaller values of 
p, as a first approximation, one may neglect the contribu- 
tion with large p and focus on the p = limit. Note that 
in the case of multiple emissions leading to the change of a 
fragmentation function, the momentum fractions carried 
out by the gluon is restricted to a small fraction of unity. 
This consideration adds additional weight to the neglect 
of these terms, especially for the case of leading hadron 
suppression calculations. We may further shift the l q , in- 



tegration as I, 



9_L 



9_L 



■l±, to represent the difference in 
the transverse momentum of the radiated gluon and the 
final quark. Incorporating the above mentioned simplifi- 
cations, expanding the transverse momentum dependent 



dY ± \k' l ± \ 2 e^-y^(p\A+(y ± /2)A+(- y± /2)\p) (81) 



d 2 y le^^{p\F+^(y ± /2)F+J-y ± /2)\p). 



In the equation above, F + ' ai represent the gluon field 
strength operators. It should be noted that only runs 
over the transverse components. In the full expression for 
Eq. , the incoming transverse momenta k 1 ^ , p 1 , appear 
only in the exponentials and may be integrated out to 
confine the transverse locations to the origin. Hence we 
define the longitudinal position dependent coefficient, 



D{y~) 



2-K l a. 



d5y- 



d 2 k_ 
! (2tt) 2 



X d 2 y\^ yi ± (82) 



x (p\F+»(yr + SyT,y\/2)F+{yT,- y y2)\p) 

y *-(p\F+^yr+5yr)F+(5yi)\p). 



2iT 2 a, 



2n 



As the integration over the transverse momentum k l ± re- 
stricts the transverse separation y 1 ^ to zero, this is no 
longer explicitly written out in the final expression. 
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This last simplification allows one to express the cross prior to emission (i.e., restricting to only the first 6- 
section for N symmetric scattering with no scattering function of Eq. ([70]) with p — > 0) as, 



mN,N 
u 0,0;m 



dyd?l± 2 4na s C F P(y) 
(2tt) 3 q ~ 



JV-1. 
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(83) 



Summing over m for fixed iV combines the gluon and 
quark color factors, as 



N\ 



N 

m=0 y ' 



C m C N-m = + Cp j 



N 



(84) 



As in the case of photon production we absorb the color, 
phase and density factors into the definition of two new 
diffusion coefficients: 



D = 



E(x L y- 



(C F + C A ) P D 

(2 P +(N*-i)y 
(c F + c A )p / dcr 



(85) 
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2p+(/V 2 - 1) 



Summing over N allows us to resum the multiple 
derivatives of the two dimensional delta function as [371 ] . 

- ( [ dy- D(y-)) N (V 2 ) N 



N=0 



(86) 



4tt / dy-D{y-)' 



Using the above identity one may easily resum Eq. 
over all scatterings. The final result is a very simple 
extension of the result of next-to- leading twist. 



B. Color factor and length enhancement for more 
complicated diagrams 

After demonstrating the outcome of a rather simple 
and leading case, the color and length factors of more 
complicated diagrams will be outlined in this subsection. 



The entire set of diagrams may be divided into six dif- 
ferent types with different color and length enhancement 
factors. Any diagram may be constructed as a combina- 
tion of these six types of diagrams. The first and simplest 
type of scattering are those which occur before the gluon 
is radiated, as shown in Fig. [4] All scatterings on the 
amplitude and complex conjugate are correlated in pairs 
and ordered in time. The general color factor for p scat- 
terings is given as, 



Tr 



i=l 



t b '«ct a ° Yl t° 
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JVc 



n 

»=i 



5 a * 



A 2 - 1 



-Tr 



[t b « 



Ct a 



(87) 



A c (iV 2 - 1)- 

As all the scatterings are correlated and ordered the 
diagram has the obvious length enhancement factor of 
{y p ) p jp\ where we have assumed that the quark starts at 
the origin and the p th scattering occurs at the location 
y p . The N c x N c dimensional operator C represents fur- 
ther scattering encountered by the quark and gluon pair 
after emission. This particular type of scattering has no 
color entanglement effects with other types of scattering. 

The next type of scattering diagrams are sketched in 
Fig. [5] and involve the interference of scattering that oc- 
curs before the emission of the gluon correlated with scat- 
tering on the quark produced after the emission. The 
color factor for this type of diagram is given as 



Tr 



wJl 



N, 



W 1 



(C F - C A /2) r 



Tr 



[t b '» 



Ct° 



n c (n? - xy 

Very similar to this type of diagram is the case where 
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FIG. 4: Simplest example of scattering prior to emission. FIG. 6: An interference effect where the scattering of the 

quark before emission is correlated with the scattering of the 
ghion after emission. 




FIG. 5: An interference effect where scattering of the quark 
before the emission is correlated with scattering of the quark 
after the emission. 



the scattering on the quark before emission are correlated 
with scattering on the gluon after emission. The color 
factor for this type of diagram, shown in Fig. [6] is given 

as, 



Tr 



n 



l[(-it<)t b °6t 



i=l 

N* - 1 



1 m 

I I fbi-idibi xb a Q Ab m b Q 



{C A /2) m Tr \t b ™t b '°C 



N C (N? - 1)' 



cf b '°. (89) 



Since these scatterings are also ordered and correlated, 
we find that the length enhancement factor in this case 
for m scattering is simply (y m + p — y p ) m /m\ where y p is 
the location where the scattering begins and y p + m is the 
location of the last scattering of this type. The color con- 
tributions from such diagrams are completely untangled 
with the case of Figs. @] and [5] and also with later scatter- 
ing on the quark and gluon after emission. The operators 

C and (j\ b ° indicate further scattering. 

The remaining three types of scattering are all in the 
final state. Two of these are rather simple and involve 
correlated ordered scattering of the final quark or the 



gluon which yields the obvious multiplicative color fac- 
tors of C$/{Nl - 1)™ or C\/{N% - l) n and have no inter- 
ference with each other yielding separate length enhance- 
ment factors of (yf — y p + m ) n /nl. As these are final state 
effects, yf is the extent of the medium and y p + m is the lo- 
cation of the last interference term between scattering in 
the initial and final states. As these color and length fac- 
tors are rather trivial no diagram corresponding to these 
is presented. 

The most complicated situation, which is ignored by 
most other formalisms is the interference between quark 
scattering in the final state and gluon scattering in the fi- 
nal state along with entangled diagrams where such scat- 
tering interchanges with ordered correlated scattering on 
the final quark and gluon. Analyzing the purely interfer- 
ing case of the type of diagrams shown in Fig. [7] with n 
scatterings, the length enhancement factor may be imme- 
diately surmised to be of the form (yf — y p+m ) n / n\ . Note 
that these contributions and those of scattering purely 
on the final quark and the gluon may become entangled. 
However at the same order of coupling, the former will 
always yield a larger length enhancement factor in the 
limit of multiple scattering. 

The color factor for diagrams of the type shown in 
Fig. [7J may be calculated as, 



C n = Tr 



f[(-it<)f° 



1 " 

— TT/ C 



6 a ' a i 

m - 1 



(90) 



There is no general formula for terms of this type. For 
the case of n = 1, we obtain, C\ = Cf ^2^( 2 \ f° r n = 2 
we obtain C*2 = 0. For n = 3 we obtain, C3 = jfj^r^fyz- 
Further scatterings cannot be done analytically and nu- 
merical methods yield that these color factors are very 
suppressed compared to all other types of scattering at 
the same order of coupling. These will henceforth be 
ignored in the remaining analysis. Diagrams with in- 
terchanging attachments where the first scattering on a 
quark in the amplitude interferes with a gluon in the com- 
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FIG. 7: An interference effect where the scattering of the the 
quark after emission is correlated with the scattering of the 
gluon after emission. 



plex conjugate followed by a scattering on the gluon in 
the amplitude interfering with a scattering on the quark 
in the complex conjugate will yield similar color factors. 
There remain the case of entangled diagrams, which may 
be shown to be suppressed as well and will also be ig- 
nored. 



VI. THE MULTIPLE SCATTERING 
RESUMMED SINGLE GLUON EMISSION 
KERNEL. 

The analysis of the previous sections will now be in- 
corporated to obtain the leading result for the all-twist 
resummcd single gluon emission kernel. At this junc- 
ture we will also attempt to connect with other leading 
formalisms for jet energy loss. In particular we will dis- 
cuss two limits of the calculation which deal with the re- 
sults of Refs. [|| (referred to as the Arnold-Moore- Yaffe 
(AMY) formalism) and those of Refs. Q (referred to as 
the Armesto-Salgado- Wiedemann (ASW) formalism). 

We begin this final section with a discussion of two 
types of diagrams which have not yet been discussed. 
Virtual contributions for which the cut line only goes 
through one side of the radiated gluon without cutting 
the radiated gluon or the final quark. These virtual cor- 
rections are important for inclusive observables and tend 
to unitarize the sing le gluon emission kernel. In single 
hadron inclusive [25( (and also for parts of multi-hadron 
inclusive calculations (36j ) these diagrams are incorpo- 
rated by insisting the product of the final fragmentation 
function and the integral over the integrated splitting 
function be subtracted from the real contribution. In 
the case of calculations of the energy lost in a single 
gluon emission or the change in the partonic distribu- 
tions, which are more exclusive i.e., depend on the mo- 
mentum carried out by both quark and gluon, these dia- 
grams play no role as one is more interested in solely the 
real contributions. 



The other type of diagrams are those which have an 
unequal number of scattering on either side of the cut in 
either real of virtual diagrams. Such diagrams unitarize 
the effect of multiple scatterings. As shown in Ref. fl2T ]. 
summing over all types of multiple scattering and inte- 
grating over the transverse momentum of a given line 
should yield a factor of order unity i.e., no length en- 
hancement. This effect is introduced into the calculation 
by seeking normalized solutions of the diffusion equation 
that results from the multiple scattering expansion. 



The limit of only gluon scattering and leading 
twist 



The is essentially the limit of large N c , where we as- 
sume that only the gluon scatters. As demonstrated by 
the color factors derived in the previous section, the dif- 
ference between a quark and a gluon scattering in the 
medium induces a factor of Ca/Cf = N^/[2(N^ — 1)] 
per scattering, which is enhanced n times for n scattering. 
Thus in the limit of multiple scattering, the naive expec- 
tion is that such diagrams dominate over those where the 
quark also scatters. This is not entirely accurate as for 
a given order the length enhancement factor for the case 
of equal scattering on the final quark and the gluon is 
length enhanced over the case of scattering solely on the 
quark or the gluon. The actual situation is intermediate 
between these two cases. In the small y limit, the scat- 
tering on the initial quark is suppressed by y. Also for a 
computation of the energy loss, the scattering on the final 
quark does not yield any change in the momentum de- 
pendent part of the expression for the energy loss. This is 
because, no factor of the transverse momentum imparted 
to the final quark appears in any of the terms of the in- 
tegrand except for the transverse momentum (5-function. 
The scattering on the final outgoing quark however does 
introduce a non-trivial change in the phase factor which 
will be considered in the next subsection. 

One may easily obtain the expressions corresponding 
to the limit of only gluon scattering by setting m = N 
in Eq. ([551) . All factors of Cp which arise from the scat- 
tering of the quark disappear. The sole factor of Cf cor- 
responds to the gluon emission vertex. As a result, the 
only type of transport coefficient that plays a role is the 
gluon transport coefficient q, i.e., the diffusion coefficient 
of Eq. (|55j) simplifies to 
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The transverse broadening of a parton which has pro- 
gressed a distance L~ , is given as 



Zi=4/ dCD(C) 



dC 



d(q. (92) 



The factor of 1/2 along with q is due to our notation of 
light cone variables where £~ = 2Q and the definition of 
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q as l\/L. Thus we obtain, 



q = 8D = 



8Tr 2 a s C A 
N 2 - 1 



px T G(x T ) 



(93) 



the transverse momentum (5-function may now be re- 
summed. To be able to resum the result, one has to 
extract the first length integration over and order the 
subsequent scatterings. As a result, we obtain, 
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N-l 



dy-D(y-) 



In order to sum over N, and indeed even in the derivation of the expression above for a fixed N, we have ignored 
terms which are suppressed by powers of N in the large N limit. While this allows us to resum the series above into 
a closed form, it has the disadvantage of not reducing to the known limit at N = 1 [TJ| ■ Summing over N we obtain 
the single gluon emission kernel in the ASW limit of multiple scattering as, 
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(95) 



The above expression represents the leading contribution 
at all twist to the double differential distribution of a 
hard quark and a gluon produced from a hard initial 
quark. While we keep referring to this as the ASW limit, 
it should be pointed out that we are evaluating this in 
the range l\ 3> k 2 ^ to l\ > kj_, whereas in the original 
ASW papers, the evaluation is carried out toward the 
lower limit. 



Guided by our assumption that l± is large, we have 
expanded the complete result of Eq. (|70p as a series in 
kj_/l± keeping terms up to l/lj_. Expanding further will 
lead to more terms. However, the Eq. (|95p has a special 
property when applied to inclusive observables such as 
energy loss. To derive the single gluon emission kernel 
from the double differential distribution above, we inte- 
grate over l q± and over the angle of gluon emission. One 
will immediately note the identity: 
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as well as the vanishing contributions, 



d %±Vl 9± 7Z r J -- T ,,..^ = °. ( 97 ) 
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The second identity is true only in the case of unpolarizcd 
observables. Integrating over l q± , we obtain the simple 
and well know result for single scattering induced emis- 
sion P3, 
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(98) 
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While the above result may come as a surprise, the 
primary reason for the cancellation of the large number 
of terms is the unitarity of the transverse momentum 
distribution due to multiple scattering. The sum over di- 
agrams, where the scatterings are not nested as in Fig. El 
restores unitarity to the distribution and thus enforces 
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the choice of a normalized Gaussian solution to the dif- 
fusion equation [Eq. ([55]) ]. In obtaining the above re- 
sult we ignored higher transverse momentum derivatives 
which contained larger powers of l\ in the denominator. 
Including these will introduce corrections to the above 
simple form. This is tantamount to restricting the cal- 
culation to the next-to-leading twist part of the all twist 
calculation. As would have been expected this yields the 
same result as the small y limit of the single scattering 
calculation [Toj . 



B. The full calculation and comparison with the 
AMY and ASW schemes 



In this subsection, we incorporate the effect of multiple 
scattering on both the quark and the gluon in the calcula- 
tion of the single gluon emission kernel. In the preceding 
section, the various types of scattering were divided into 
six categories. Of these, two types which involve the in- 
terference between the scattering of the gluon and the 
quark after emission have very small color factors and 
do not possess larger length enhancements compared to 
other diagrams. These will be ignored in this effort. This 
leaves only four different types of diagrams that need to 
be considered: scattering on the quark before emission, 
interference between scattering of the quark or gluon af- 
ter emission with scattering before emission and separate 
scattering of the produced quark and gluon. Diagrams 
which entangle these contributions tend to be suppressed 
in length factors. Thus these four contributions can be 
considered one after the other separately. 

In order to incorporate all these effects. We divide 
the length of the medium traversed from y^ to L~ into 
four regions. The first part from y§ to yq within which 
only the initial quark scatters and the gluon emission has 
not occurred in either amplitude or complex conjugate 
(FigU). From yc > yQ to yE where there is interfer- 
ence between the scattering of the quark in the ampli- 
tude (complex conjugate) with scattering of the quark 
or gluon in the complex conjugate (amplitude) [Fig.[5]or 
Fig. E] . These are referred to as the cross terms (hence 
indicated as yc) and continue till the later gluon emis- 
sion at yE- We also keep track of the location (c which 
is the first scattering of the gluon in the complex con- 
jugate (amplitude). The last region is from yi > yE to 
L~ which represents the part where the produced quark 
and gluon scatter independently. We also keep track of 
Cf ^ Ve which is the first scattering of the gluon in this 
regime. Note that Cc < Cr < L~ . The gluon emission 
kernel will depend on all these locations. 

Unlike in the preceding subsection, where only one dif- 
fusion coefficient was required as there was only one type 
of scattering (that on the final state gluons), there are 
four different types of scattering in this subsection. As 
a result, four different types of diffusion coefficients will 
be required which differ in their overall color factor only. 
We denominate these as Dqg (or simply Dq) where both 



scatterings in the amplitude and complex-conjugate oc- 
cur on the gluon, 
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The standard transport coefficient in jet-quenching q = 
8Dgg- The case where one scattering is on the gluon 
and one on the quark is denoted as Dqq or Dqq inter- 
changeably, and given as 



D 



GQ 



{C A /2)pD 
2p+(A c 2 - 1)' 



(100) 



The case where both scatterings are on the quark and 
one occurs prior to gluon emission and one occurs after 
will be designated as Dqq and has the expression, 



D 



QQ 



(C F - C A /2)pD 
2p+(7V c 2 -l) ' 



(101) 



It may come as a surprise that this diffusion coefficient 
is actually negative. Multiple such insertions of inter- 
fering terms where the quark in the amplitude scatters 
before the emission and the quark in the complex conju- 
gate scatters after the emission lead to a shrinking of the 
transverse momentum between the radiated gluon and 
the final outgoing quark. The final diffusion coefficient is 
the case where the the quark scatters in both amplitude 
and complex conjugate and both happen either before or 
after the gluon emission. The diffusion coefficient in this 
case is designated as Dqq, or simply Dq, and given as 



Dqq = 



(C F ) P D 
2p+(N 2 - 1)' 



(102) 



The derivation of the full double differential distribu- 
tion of the cross section as a function of the transverse 
momentum of the quark and the gluon is rather compli- 
cated and will be presented elsewhere. In this article, we 
will simply derive an expression for the energy loss dis- 
tribution. From the analysis of the preceding sub-section 
on the ASW limit we note that in the differential distri- 
bution, all terms which involve extra derivatives of the 
delta function vanish in the integrated expression. Also 
terms which involve one integral to be kept separate from 
the rest due to a phase factor depending on this location 
will primarily be involved in the higher twist contribu- 
tions. Thus terms with more than one phase factor are 
suppressed by larger powers of l\. Thus to obtain the 
leading contribution only the terms with one phase fac- 
tor need to be retained. 

In the interest of simplicity we will not derive the extra 
contribution from diagrams where the quark scatters be- 
fore the emission without interfering with any other type 
of scattering. These will be included separately at the 
end of the section. These diagrams were shown to have 
smaller length enhancement factors. They proportional 
to y 2 and are thus very suppressed in the small y limit. 
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As a result, they play little role in the computation of 
the medium modified fragmentation function [TTj where 
the range of the gluon momentum fraction is limited by 
< y < 1 — z. The fraction z is the momentum fraction 
of the detected hadron with respect to the initial hard 
quark's momentum and is thus also the minimum mo- 
mentum fraction that must be carried by the final quark 
to be able to fragment and produce such a hadron. These 
diagrams also do not contain a vacuum like contribution 
in either amplitude or complex conjugate and arise as a 
purely medium dependent piece. As a result, the neglect 
of these terms focuses this part of the calculation to man- 
ifestly involve the vacuum contribution in either the am- 
plitude or complex conjugate. As a result, this part of the 
calculation is somewhat orthogonal to the calculations in 
the formalism of Arnold, Moore and Yaffc (AMY) where 



none of the terms contain a vacuum contribution. This 
part of the calculation, however has a considerable over- 
lap with the calculation of Armesto Salgado and Wiede- 
mann (ASW) except that in the ASW case there is no 
Taylor expansion in the exchanged momentum i.e., there 
is no explicit assumption lj_ kj_. 

Invoking the above mentioned simplifications, we write 
the expression for g^s+mm " where q = r+s. This 
represents the operator expression for the case where the 
quark emits a gluon at the origin in the amplitude and 
after q scatterings in the complex conjugate. The q scat- 
terings interfere with r scatterings on the produced quark 
and s scatterings on the gluon in the amplitude. After 
these interfering terms, the produced quark and gluon 
scatter m and n times in both amplitude and complex 
conjugate. The full expression is written as, 
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As mentioned above, we will consider the case where the initial state scattering on the quark is ignored thus we set 
= 0. The least power suppressed contribution in this limit is the term containing the e lp XLy ° = 1 factor. 
The second line in Eq. (|103l) may be simplified as, 
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In the equation above, we have only retained only terms which are suppressed by at most l\. This is all that is 
required for the case of terms containing the e %p XLV ° = 1 factor. Terms containing two separate phases will require 
the products (£| =1 fcj_) 2 (£™ i ^l) 2 which appear in the above expansion suppressed by While it is not difficult 
to include these terms and these will be included in future numerical calculations, we refrain from including them 
here as we illustrate the leading behavior of the multiple scattering gluon radiation spectrum. 

Expanding the 5-function and the (non-phase) momentum dependent part as a series in transverse momentum, we 
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obtain, 
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In the equation above, the g integrals run from yo = 
to £j or y j which ever comes first. The q integrals are all 
ordered and r of them attach to a quark on both sides 
of the cut. For terms without the exp[—ip + XLy E ], this 
yields the length enhancement factor of, 
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For the term containing the y E phase factor, if y F scat- 
ters on a quark on both sides, then we take the r th inte- 
gration out from the product. If y E scatters on a gluon 
then we take the s th integration out from the ordered 
product, yielding the factor, 
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As a result, including all the length dependent terms, 
associated with the y E phase factor, we get, 
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When combined with the transverse momentum deriva- 
tives one notes that there is no contribution from the 
terms containing r — 1 (as these will always contain an ex- 
tra Vf which is not being summed) i.e., where y E repre- 
sents the scattering on a quark. Summing over r, s, m, n, 



and integrating over l q± , the only surviving term is, 
0(l)=/^«.CP(y) 



x / dy E AyD GQ {y^)e-^ (109) 
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Note that this term is explicitly y dependent and thus 
subleading as y —> 0. 

Similarly for terms with the e~ lp Xh ^ 1 phase factor, 
the length dependent terms, after summation over all 
insertions and integration over l q±1 yields, 
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The two equations above, represent terms which are 
not squares of amplitudes but rather interference terms. 
Thus we also need to include the complex conjugates. 
The result of summing the two equations above with their 
complex conjugates may be easily obtained by converting 



XlV e to 2cos(p + XLy E ) in Eq. (|109j) and e 



-ip + x L C 

to 2cos(p+;elC~) in Eq fTTU]) . 

In the results above, we have only obtained terms 
which contain a cosine of the usual phase argument 
{xlp + Q~) but not the terms proportional to unity. The 
combination of these yields the LPM effect. To obtain 
the terms which don't contain the argument xlP + (~ , we 
have to consider terms with at least two of the sums re- 
stricted, i.e., terms where r = 0,s = or terms with 
m = 0,?i = 0. Considering these latter contributions 
changes the overall phase factor, e.g., when r = 0, s = 
one will have to replace with yj and £^ with £j~ in 
Eq. (|103[) . While this limit provides the zeroth order 
contributions to the various diffusion terms that go with 
the cos(p + xlC _ ), it also introduces new contributions, 
i.e., the terms which don't depend on the cosine, which 
have to be included separately. This we carry out in the 
following. 

The reader will note that all the sums over the number 
of scatterings start from a minimum of 1, e.g., in the first 
line of Eq. (|103[) . This requires that all types of scatter- 
ings considered have to have taken place at least once. 
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We also need to include the cases where a particular type 
of scattering did not take place. The simplest situation is 
the case where m = n = r = s = 0. In this case there will 
not be any of the phase factors or in-medium matrix ele- 
ments and an inspection of Eq. (| 103|) will reveal that this 
yields simply the pure vacuum contribution of Eq. (|17[) . 
The second case is where we set r, s the cross scattering 
to 0. Beyond all other trivial changes, this modifies the 
phase factor to yield, 
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We ignore the contributions that contain yj as there 
is no corresponding p± that appears in the amplitude, 
thus these will vanish when integrated over l q± . Also 
all terms which contain the phase factor exp (±ip + xlCT) 
give leading contributions when yj > C7 • This allows the 
integrals over quark scatterings to be lumped together 
and included in the broadening of the two dimensional 
(^-function over l q± . Thus the leading contributions (i.e., 
the least suppressed in Q 2 ) which contain at most one 
location Q are given as, 

(l - e -* p+XL ^ - 6(yJ - Ci)er lp+XL ^^j 
X (l - e ip+XL( ' - 9{yJ - CT)e ip+XLC ' 

= 2 + 3%7 - Cj") - 2 cosO+^Cj") 
- 2%7-C7)cos(p+z L C7)- (H2) 

Along with these one will have to consider the terms 
which depend on the location yj , i.e., 
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There is one extra factor of 9(jjj — Cj ) without a cosine 
in Eq. (|112[) and one less factor if 0(C7 — 2/7) (without 
a cosine) in Eq. Q113p . Since these terms have no extra 
dependence on any position they are essentially equal, 
i.e., these are symmetric across the line yj = As a 
result, the term containing the extra factor or 9(yJ — Qj) 
can be converted to a term containing an 9(yJ — (y) and 
included with Eq. (|113[) . Such a conversion is required in 
the l± 3> k 1 ^ limit, as any term which cannot be cast in 
the final form, 
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up to multiplicative factors, diverges as at small l±. 
Such a divergence calls into question the approximation 
of Taylor expanding in the softer momenta fcj_. Such 



a term also dominates over the vacuum radiation term 
at small l± and thus will invalidate the use of an in- 
medium DGLAP evolution equation to incorporate the 
effect of multiple radiations. All terms, at this order of 
approximation, are finite at l±_ — > due to the destructive 
interference of the LPM effect. 

We now consolidate terms which contain one more in- 
tegral over light-cone position C^J than transverse deriva- 
tive on the S function, in Eq. (|105|) . These are of two 
types, those which remain as n — > and those which 
vanish in this limit. Terms of the first type may be ex- 
pressed as, 
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Summing over m, n reduces the second line of the above 
equation to the broadened Gaussian, 
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The terms in Eq. (|105[) which only exist for n > 1 
and contain one more integral over the light-cone location 
than transverse derivative on the d function can be 
expressed as, 
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To sum over m, n one of the n integrals have to be ex- 
tracted from the sum. Denoting this location as y^ , the 
result of the sum over m, n for the second line of the 
above equation yields, 
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Integrating the above term over l q± yields zero due to 
Eq. dST]). 

Combining these contributions, yields the leading 
medium dependent contribution to the gluon emission 
spectrum from multiple scattering in the small y limit 

as, 
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Note that this is identical to the contribution in the pre- 
vious subsection where we insisted that only the gluon 
scatters. 

The next-to-leading contribution in the small y limit 
may be derived similarly from the contributions in 
Eq. ([109[) . including in addition the contributions which 
do not contain an exponential or a cosine factor, as 
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r)] 



We now add the contributions where the initial quark 
scatters independently, i.e., where none of these scatter- 
ings interfere with scatterings on the radiated gluon or 
on the final quark. These diagrams contribute terms sup- 
pressed by one extra power of l\ only in the limit where 
r, s = 0. These may be written down in analogy with 
Eq. (|119[) by replacing k l ± — > yp* ± and Dq with Dq. As 
a result, this y suppressed contribution is given as, 
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While the expressions above may be used in an en- 
ergy loss calculation to determine the medium modified 
fragmentation function, it cannot be use for the Monte- 
Carlo simulation of jet propagation in medium. That 
will require the full double differential expressions in lj_ 
and l q± . While these are not presented here, they may 
be straightforwardly derived from Eq. (|103j) and will be 
discussed elsewhere. 

Though not discussed in this article, the calculation 
of the medium modified fragmentation function or in the 
Monte-Carlo simulation of a hard jet, one requires proba- 
bility decreasing contributions which counter the effect of 
the probability increasing expressions such as Eq. (|1 191) . 
In the case of the medium modified fragmentation func- 
tions, these are essentially introduced as unitarity (or 
virtual) corrections. In the case of jet simulations, one 
will traditionally derive a Sudakov factor [38{ from the 
real splitting correction. 



VII. DISCUSSIONS, CONCLUSIONS AND 
OUTLOOK 

In this article, we have extended the higher twist for- 
malism by deriving the leading length-enhanced contri- 
bution to the single gluon splitting kernel at all-twist. 
Physically this includes the process of gluon emission 
from a hard quark which undergoes multiple scattering 
within the formation time of the produced collinear ra- 
diation. 



We considered the propagation of a hard quark pro- 
duced in DIS on a large nucleus. The hard virtual photon 
strikes a hard incoming quark in a nucleon which then 
proceeds to burrow through the nucleus in the direction 
of the hard photon. The hard quark is off-shell at the 
point of production and in the absence of the remaining 
nucleons would proceed to radiate multiple gluons prior 
to hadronization. We focused on one such hard emis- 
sion l± 3> Aqcd which occurred inside the nucleus. The 
presence of the medium modifies the cross section to ra- 
diate the hard gluon due to multiple scattering of the 
propagating quark and the radiated gluon. 

The main focus in the current effort has been on 
the calculation of the triple differential cross section 
to produced a single gluon with transverse momentum 
momentum fraction y and a final outgoing quark 
with transverse momentum l q± (and momentum fraction 
1 — y). The series of diagrams which contain the leading 
length enhanced corrections to this process were identi- 
fied. These diagrams were then evaluated in the collinear 
and hard approximation q~ l±_ 3> fc_i_ where the first 
inequality suggests that the jet energy is logarithmically 
larger than the transverse momentum of the radiation, 
which by the second inequality is power enhanced com- 
pared to the transverse momentum imparted from the 
medium. 

Incorporating the various simplifications that result 
from the limits mentioned above, we evaluated the lead- 
ing terms in the triple differential cross section. Besides 
the effect of multiple scattering on vacuum like splitting 
diagrams, we also obtained new contributions which rep- 
resent the physical effects of a hard incoming gluon hit- 
ting a near on shell quark and the forcing it off shell and 
inducing a radiation. This interferes with the vacuum like 
emission as well as with other processes where on-shell 
quark splits, radiating a gluon and either the out-going 
quark or gluon goes off-shell. The off-shell parton then 
quickly loses its excess virtuality by encountering another 
collision with an incoming hard gluon. 

After a discussion of length enhancements and color 
factors of different types of diagrams, we focused on di- 
agrams which yield the leading power corrections to the 
single gluon cross section. Among these, the leading 
corrections in the gluon momentum fraction y with the 
largest color factors arose from processes where the fi- 
nal produced quark or gluon scattered independently or 
where there was interference between the scattering on 
the quark prior to emission and that on the quark or 
gluon after emission. 

Following this, the double differential cross section to 
radiate a hard gluon was evaluated. This was evaluated 
in two separate cases: the case of only gluon scattering 
and the case which included all types of scattering. The 
leading corrections in a power expansion involving the 
hard scale were then isolated. In these diagrams, the in- 
sistence that the transverse momentum of the radiated 
gluon is large in comparison with the transverse momen- 
tum of the incoming gluons from the medium make this 
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calculation somewhat orthogonal to that of AMY. Inte- 
grating the transverse momentum of the final produced 
quark reduced these contributions to those obtained from 
a single scattering analysis. This implies that the lead- 
ing effect of multiple scattering is to broaden the trans- 
verse momentum distributions of the produced quark and 
gluon leaving the energy loss kernel unchanged. This jus- 
tifies a similar assumption make in Ref. (39j , where it was 
assumed that the broadening of a radiated gluon could 
be considered independently from the emission cross sec- 
tion. 

In this paper we have restricted the calculation to only 
the analytical component. The numerical evaluation of 
the medium modified fragmentation function followed by 
the phenomenological implementation in DIS on a large 
nucleus and in jets in high energy heavy ion collisions 
will be carried out in a separate calculation. These will 
include the effect of multiple such stimulated emissions. 

Along with the various terms outlined in the preced- 
ing section, we will also include a series of other terms 
which were not explicitly evaluated in this effort. Princi- 
ple among these are the contributions suppressed by Q A 
(or rather and greater. The neglect of such contribu- 
tions is justified at large Q 2 . For a Q 2 ~ 3 — 5GcV 2 , as 
is the case for the HERMES experiments, and for large 
enough nuclei, such power corrections may also have to 
be included. 

The DGLAP equation in vacuum, is a coupled equa- 
tion which connects the evolution of the quark fragmen- 
tation function with that of the gluon. This will also be 
the case for energy loss in a medium. Thus a complete 
in-medium evolution will also require a calculation of the 
in-medium splitting function for a gluon to two gluons or 
a quark anti-quark pair. The derivation of such a split- 
ting function will proceed in the same manner as in this 
article and will be presented in a future effort. 
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IX. APPENDIX A: A SCALING OF THE A+ 
FIELD 

We use the linear response formula to ascertain the 
power counting of the A + field. Suppressing the color 
superscripts we obtain, 



A»(x) = / d A yV^(x-y)J v (y). 



(122) 



In the equation above, D is the gluon propagator and at 
leading order in the covariant gauge is given as, 



d A k 



(2tt)< 



k 2 



le 



(123) 



In Eq. (|122p . J v {y) = -ip(y)j u ^{y) is the current of par- 
tons in the target which generates the gluon field. The 
fcrmionic operator may be decomposed as, 



ip(y) = 



dp + d 2 p± 



(2tt)\/ P + 



vp-y 



(124) 



The scaling of the fermionic operator depends on the 
range of momentum which are selected from the in-state 
by the annihilation operator. Note that this influences 
both the scaling of the annihilation operator a p as well 
as the bispinor u(p). The power counting of the annihi- 
lation operator may be surmised from the standard anti- 
commutation relation, 



{a r v ,a s ^}={2-Kf5 i {p-0)S 



(125) 



and the power counting of the bispinor from the normal- 
ization condition, 



7 p + + 7 + p 



7-L • PX- 



(126) 



Substituting the equation for the current in Eq. 11231 and 
integrating out y, we obtain, 



A' 



d 3 pd 3 q 



-i(p-q)-x 



(P - qf 



-a\a p u(q)~f + u(q\l27) 



If the incoming and out going momenta p and q scale 
as collinear momenta in the (+)-direction, i.e., p ~ 
Q(1,A 2 ,A), then we get, S 3 (p-p l ) ~ [A 2 Q 3 ] -1 , as one 
of the momenta will involve the large (-f)-component 
and the remaining are the small transverse components. 
Thus the annihilation (and creation) operator scales as 
\~ 1 Q~ 3 / 2 . Also in the spin sum p 1 ~ Q and thus 
u(p) ~ u(q) ~ Q 1 / 2 ; one can check that the 7 + projects 
out the large components in u and u in the expression 
u{q)^ + u{p). We also institute the Glauber condition that 



a- 



V 



A Q and p±_ — q±_ ~ \Q. 



Using these scaling relations we correctly find that the 
bispinor scales as AQ 3 / 2 . However, to obtain the correct 
scaling of the gauge field A + one needs to institute the 
approximation that q + = p + + k + where k + ~ A 2 Q. 
This condition is introduced by insisting that the (+) 
momentum of the incoming and outgoing state, which 
control the scaling of a' q and a p , are separated by k + ~ 
A 2 Q. This is used to shift the dq + — S- dk + and as a result 
we obtain the scaling of the A + field as A 2 Q. Following 
a similar derivation, with the replacement j + — > 7 we 
obtain the scaling of the A 1 - <~ A 3 Q. 
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